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ABSTRACT 


A  technique  for  using  a  sequence  of  two-pulse  bursts  for  obtaining 
estimates  of  mean  wake  velocity  and  wake  velocity  width  has  been  proposed 
by  W.  D.  Rummler  of  Bell  Telephone  Laboratories.  The  derivation  of 
this  technique  was  based  on  heuristic  reasoning.  The  present  note  derives 
the  same  estimates  from  the  theory  of  maximum  likelihood  estimation. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 


In  a  recent  series  of  memoranda,  W.  D.  Rummler  of  the  Bell 
Telephone  Laboratories  has  proposed  simple  estimates  of  mean  wake  velocity 
and  wake  velocity  width.  These  estimates  were  obtained  on  the  basis  of 
purely  heuristic  arguments.  The  purpose  of  the  present  note  is  to  show 
that  these  same  estimates  can  be  derived  from  fundamentals  using  the  method 
of  maximum  likelihood. 

The  transmitted  signal  s(t),  to  be  used  for  obtaining  these  estimates 
is  a  sequence  of  two-pulse  bursts  such  as  shown  in  Figure  1.  The  pulse 
width  A  is  chosen  to  yield  the  desired  range  resolution  and  the  interpulse 
spacing  is  chosen  so  that  the  wake  returns  from  the  two  pulses  of  the  pair 
do  not  overlap.  The  time  between  bursts  is  made  large  enough  so  that  the 
returns  from  a  range  cell  located  a  fixed  distance  behind  the  re-entry  vehicle 
are  statistically  independent.  This  condition  of  statistical  independence  can 
also  be  achieved  by  suitably  jumping  the  carrier  frequency  between  successive 
two-pulse  bursts. 

The  signal  +  r(t),  received  when  a  single  burst  is  transmitted  will  be 
taken  to  be 


r(t)  =  c(t)  +  n(t) 


(i) 


where  n(t)  is  white,  Gaussian  receiver  noise  and 


c(t)  =  JJa(T,f)s(t-r)ej2lTftdTdf  (2) 

is  the  uncorrupted  wake  return.  The  function  a(r,  f)  is  the  complex  amplitude 
density  of  that  portion  of  the  wake  return  having  delay  j  relative  to  the  re¬ 
entry  vehicle  and  Doppler  shift  f.  This  function  is  random  with  the  properties. 

E[a(r,f)]  =  0 

E[a  (t,  f)a(T  f 7 )]  =  0 

E[a(T,f)a*(T',f')]  =  a  (t,  f )5  (r-T '  )6  (£-£' ' )  (3) 


+A11  time  functions  will  be  understood  to  be  the  complex  envelopes  of  their 
corresponding  real,  r.  f.  signals. 
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Fig.  1.  The  two-pulse  bursts. 
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where  a(r,f)  is  the  average  scattering  cross-section  density  of  that  part  of 
the  wake  located  at  delay  t  and  Doppler  shift  f.  The  basic  problem  at  hand 
is  to  estimate  the  nf-centroidM  and  nf-width,f  of  o(T,f)  for  various  values  of  t. 

The  conditions  (3)  correspond  to  assuming  that  the  returns  from 
different  wake  scatterers  are  statistically  independent  and  that  the  phases  of 
the  returns  are  uniformly  distributed.  It  will  further  be  assumed  that  the 
total  number  of  scatterers  is  large  enough  so  that  some  form  of  the  Central 
Limit  Theorem  holds  and  c(t)  can  be  modelled  as  a  Gaussian  process.  This 
being  the  case,  c(t)  can  be  characterized  completely  by  the  three  functions. 

E[c(t)]  =  0 
E[c(t)c(t')]  =  0 


E[c(t)c*(t')]  =  JT<7(t,  f)s(t-T)s*(t,-T)e^2lT^t”t  Wdf  (4) 

In  order  to  facilitate  the  forthcoming  analysis,  it  will  be  assumed 
that  the  radar  receiver  operates  in  a  range- sampled  mode.  This  means 
that  the  received  signal  is  passed  through  an  i.  f.  filter  whose  bandwidth  is 
roughly  equal  to  the  signal  bandwidth  and  then  sampled  once  per  range  re¬ 
solution  cell  over  the  range  interval  of  interest;  in  the  present  case  this  is 
the  range  interval  over  which  the  wake  return  is  to  be  studied.  The  result 
of  this  process  is  that  each  transmitted  pulse  pair  produces  a  pair  of  complex 
samples  for  each  range  gate,  r,(i)  =  r(r  +iA),  r?(i)  s  r (t  +  T+iA)  where  t 
denotes  the  range  delay  at  which  the  first  wake  return  of  interest  occurs  and 
i  denotes  the  number  of  the  particular  range  gate  under  study.  These  complex 
samples  are  zero-mean  Gaussian  random  variables  whose  covariances  are 
given  by 

E[r1(i)r1(i)]  =  E[r2  (i)r2  (i)]  =  Efr^i):^  (i)]  =  0 

E[r1(i)r2*(i)]  =  A2a  Jct  (To+iA,  f)e"Ji'7Tf Tdf  (5) 

E[r  j(i)r  j*(i)]  =  E[r2  (i)r2*  (i)] 

=  A2a Jct (ro+iA,  f)df  +  a2 
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where  a  denotes  the  variance  of  the  range- sampled  receiver  front-end  noise. 
Equations  (5)  follow  directly  from  equations  (4)  when  it  is  assumed  that  the 
function  a  (t,  f)  can  be  regarded  a  s  constant  over  a  delay  interval  of  length  A 
so  that  the  r- integration  can  be  performed  explicitly  and  it  is  assumed  that 
a  (T,  f)  is  only  non- zero  in  a  delay  interval  of  length  less  than  T.  This  last 
assumption  is  equivalent  to  the  earlier  assumption  that  the  interpulse  spacing 
of  the  two-pulse  burst  is  large  enough  to  straddle  the  wake. 

The  pair  of  random  variables  r^(i),  r^(i)  is  Gaussianly  distributed 
with  a  density  function  given  by 


In  pfr^i),  r,(i)]  £  In  p[r(i)]  =  -  £*.(i)Aj"  £  (i)  -  In  [ti^Ia.I] 


where 


r  (i) 


’rj(i) 

r2(i) 


a  nd 


A.  = 

1 


3^+a2 

B.V 


B-V 

i  i 

n  2  ,  2 

+a 


P.2=A2Afa(T  +iA  ,  f)df 

1  w  O 


(6) 

(7) 


(8) 


g2p.  =  A2A  Pct(t  +iA,  f)e”^2lTfTdf  (9) 

11  0  o 

The  assumption  that  the  returns  from  successive  two-pulse  bursts 
are  statistically  independent  now  can  be  used  to  write  down  the  joint  prob¬ 
ability  density  function  for  the  returns  from  N  successive  bursts  as  follows, 

N  ,  2 

Inptr^i),  •  •  -£N(i)]  =  -£  (i)Ai“  £k(i)  +  ln[ir  |A.|l  (10) 

k  =  l 

The  maximum  likelihood  strategy  for  forming  estimates  of  the  mean 
wake  velocity  and  wake  velocity  width  in  each  of  the  range  cells  of  interest 
is  to  maximize  the  function  given  by  Equation  (10)  with  respect  to  these 
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parameters  and  with  respect  to  any  other  unknown  parameters.  To  accom¬ 
plish  this  end,  it  is  first  necessary  to  give  precise  definitions  of  these 
quantities.  There  will  be  done  by  assuming  that  the  function  o^r^+iA,  f)  has 
a  precisely  known  shape  but  that  its  scale,  center  and  width  are  unknown. 
Stated  in  symbols,  it  is  assumed  that 


^  f  7" 

a(r+iA,  f)  =  -ia^(i,  -zr2-) 


W. 

1 


w. 

1 


(11) 


where  the  function  a  (i,  x)  is  known  for  each  i  but  a.,  f.andw.  are  unknown 

o  111 

parameters.  It  is  convenient,  but  not  necessary  for  most  of  the  following 
analysis  to  further  assume  that  aQ(i,  f)  has  been  normalized  so  that 


fao(i,f)df  =  1 

fa  (i,  f)fdf  =  0 
J  o 

JaQ(i,  fjf2 df  =  1  (12) 

The  normalizations  defined  by  Equation  (12)  are  equivalent  to  defining  ci^, 
w^,  and  7\  in  the  manner  usually  used  in  radar  estimation  problems;  namely, 

a.  =  |  a(r  +iA,  f)df 
i  J  o 

al  =  J  a(To+iA,  f)fdf 

a.w.2  =  f  a(r  +iA,  f)(f-T.)2df  (13) 

11  w  O  1 

The  wake  parameter  estimation  problem  now  has  been  reduced  to  a 

purely  mathematical  problem,  that  of  maximizing  the  function  defined  by 

—  2  2 

Equation  (10)  with  respect  to  the  parameters  o\,  f  w„  Since  =  A  Aa^, 

maximizing  this  function  with  respect  to  a.  is  the  same  as  maximizing  it 
2  1 

with  respect  to  g.  .  This  will  be  done  in  the  following  analysis 
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Denoting  times  the  left-hand  side  of  Equation  (1)  by  Land  referring 
to  Equation  (8)  allows  Equation  (10)  to  be  written  in  the  form, 


L  =  -2y. 


(p.24a2)  -  Ree^p.x.* 

~T2 — 2  2  4, - 72 - 

(9i  -kx  )  -  p.  |p.| 


,  2r  2  2.2.4,  ,2-, 

lmr  [(g.  +a  )  -p.  IpJ  ] 


(14) 


where 


k=l  ^  '  rlk^l)'  +  '  r2k^^  1 

(15) 

,  N 

Nk-1  rlk^l)r2k"(l) 

(16) 

The  parameters  L  and  w.  enter  Equation  (14)  only  through  the  com¬ 
plex  quantity  p^.  It  is  easy  to  show,  using  Equations  (9)  and  (11),  that  any 
value  of  inside  the  unit  circle  in  the  complex  plane  can  be  achieved  by 

a  suitable  choice  of  f.  and  w.  and  conversely,  that  any  choice  of  f.  and  w. 

ii  ii 

results  in  a  value  of  p^  inside  the  unit  circle.  Thus,  maximizing  L  with 
respect  to  f  .  and  w  is  equivalent  to  maximizing  it  with  respect  to  p_^  with 
the  constraint  that  p.  lie  inside  the  unit  circle. 

The  maximum  of  L  with  respect  to  argp^  obviously  occurs  when 
argp^  =  argp^  where 


argp.  =  argx. 


(17) 


and  is  given  by 


maxL 

argp. 


=  -2y.. 


(S  .2+a2)  -  |  B  .2p  .x.  | 

i  '  i  K  i  i 

77“2  7~Z  ~F.  72 

(Bi  +a  )  -  Si  I  p  •  I 


,  2r  2  2.2  .4,  1 2, 

In  u  [(g.  +a  )  -Si  \^\  ] 


(18) 


The  maximization  of  Equation  (18)  with  respect  to  jL  and  |p  J  is  facilitated 
by  the  change  of  variables, 


f.  =  B.2(l  +  Ip.l  )  +a2  (19) 

gi  =  Si2(l  -  I  Pi|  )  +  a2 
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in  terms  of  which  Equation  (18)  assumes  the  form 


maxL 

argp.  '  "  yi 


g; 


1+  |x.| 

-  In  irg.  -  y.  — 7 -  -  In  irf. 

°i  i  f.  i 

i 


(20) 


Differentiating  Equation  (20)  with  respect  to  f.  and  g.  and  setting  the  results 
equal  to  zero  yields  the  following  equations  for  f  ^  and  g^  the  maximizing 
values  of  f.  and  g. 

l  l 

1  -  |x  1 

vi  —n - gi =  0 

g; 


1  +  |x.| 

y.  - -  f.  =  0 

yi  j  2  i 

i 

therefore, 

g.  =  y.(l  -  |x.|  ) 

=  Yj(l  +  UJ) 

and  the  corresponding  estimates  of  0^  and  |p^|  are  given  by 


(21) 


(22) 


*2  2 

p •  =  y.  -  a 

1  3 1 


Pi'  = 


nW 


yra 


(23) 


The  estimates  given  by  Equation  (23)  are  correct  (i.  e.  the  values  of 

Z  2 

3^  and  |p^|  that  maximize  L  subject  to  constraints  3^  >  0,  |pj|<  1),  if,  and 

only  if,  §i  2  >  0  and  |  p^|  <1.  If  either  of  these  conditions  is  vi  olat  ed,  then 

the  required  maximum  of  L  occurs  for  degenerate  values  of  the  parameters 

such  as  0.  =0  and/or  I  p .  |  =  L  It  can  be  shown  that  as  long  as  the  integrated 

1  21  2 

wake-to- noise  ratio  NB .  /a  is  sufficiently  greater  than  unity,  this  situation 
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will  only  occur  with  very  small  probability.  In  the  sequel  it  will  be  assumed 
that  the  wake- to- noise  ratio  is  large  enough  so  that  the  occurrence  of  the 
above  mentioned  degeneracies  can  be  ignored. 

The  results  of  the  above  analysis  now  can  be  summarized  by  com¬ 
bining  Equations  (17)  and  (23)  with  Equations  (9)  and  (11)  yielding  the  results 

2 

that  the  estimates  of  3-  and  I*,  and  w.  are  given  by  the  equations, 


arg 


2  2  2 

3  *  =  y.  -  a 

i  7 1 

.a. 

f-f . 


J-  [a  <i,^-i)e-i2'r£Tdf 

/  .  o  O  W  . 


W  . 

1 


w  . 

1 


= argx. 


i  1  f*  /•  \  -j2irfT.fi  yi^Xi' 

I  wT  K(l’^~)e  df>  = - ; 


Yi-  a 


(24) 

(25) 

(26) 


Equation  (24)  is  an  explicit  expression  for  3-  ,  whereas,  Equations  (25) 

A  1 

and  (26)  only  define  f^  and  w  .  implicitly.  These  last  two  equations  can  be 
further  simplified  by  making  the  change  of  variables 


f-f. 

i 

X  =  -7T- 


w  . 

1 


resulting  in  the  equations 

-j2irf.T 


arg 


1  K(i>x)e 


-  j2rrw .  Tx 


dx 


=  arg 


(25;) 


-  j 2ttw .  Tx  y.|x.| 

|JaQ(i,  x)e  1  dxj  ^ 


y.  -  a 


(26y) 


These  equations  still  define  the  exact  maximum  likelihood  estimates  f^  and 
w  ^  in  an  implicit  manner;  however,  they  suggest  an  approximation  that  will 
yield  explicit  expressions  for  these  estimates.  If  the  estimate  w^  is 
sufficiently  small  compared  to  (w^T  <<1),  then  the  exponential  appearing 
on  the  left-hand  side  of  Equations  (25*)  and  (26;)  can  be  approximated  by  the 
first  three  terms  of  its  Taylor  series. 
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e 


(27) 


-j2irw.Tx  ,  7  7 

a?  1  -  j2-irw  ^Tx  -  ^(2-rTW  ^T)  x 

Substituting  Equation  (27)  into  Equations  (251)  and  (26')  and  making  use  of  the 
normalizations  defined  by  Equation  (12)  yields  the  result 


f  . 

l 


1 

27 T 


arg  x 


(25") 


__1 _ 

2(ttT)2 


yra 


(26") 


Note  that  these  approximate  expressions  for  the  maximum  likelihood  estimates 

A 

f  ^  and  do  not  require  explicit  knowledge  of  the  shape  of  the  wake  dis¬ 

tribution  cr  (i,  f)  as  do  the  exact  Equations  (251)  and  (261).  Equations  (25n) 
and  (26n)  define  exactly  the  same  estimates  that  Rummler  has  proposed  on 
an  intuitive  basis.  They  are  now  seen  to  be  approximations  to  the  maximum 
likelihood  estimates  of  f  ^  and  w^  valid  as  long  as  the  estimated  wake  velocity 
width  &  ^  is  sufficiently  small  compared  to  -7^,  the  unambiguous  frequency 
interval  defined  by  the  two-pulse  burst. 

This  concludes  the  derivation  of  Rummler's  wake  velocity  estimates 
from  basic  statistical  principles.  For  a  discussion  of  the  bias  and  accuracy 
of  these  estimates,  the  reader  is  referred  to  the  references  cited  below. 
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